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Abstract. In this article, we consider Leray solutions of the Navier-Stokes equations in the exterior of 
one obstacle in 3D and we study the asymptotic behavior of these solutions when the obstacle shrinks 
to a curve or to a surface. In particular, we will prove that a curve has no effect on the motion of a 
viscous fluid, so it is an erasable singularity for these equations. 



1. Introduction 

The present article studies the influence of a thin obstacle on the motion of a three-dimensional 
incompressible viscous flow. 

Such a subject present a large literature on experiments and simulations (see the references in [7]), 
but concerning a PDE approach, we have waited 2006 with the first work of Iftimie, Lopes and Lopes 
in [6]. Therein, the authors treated the case of an obstacle which shrinks to a point in a 2D viscous 
fluid. They proved that we do not see the effect of the point on the motion of the flow, except if we 
assume a non-zero circulation for the initial velocity. Later, the present author proved in [7J that a 
curve in 2D is not erasable, as in the previous case. Finally, Iftimie and Kelliher showed in |5j that a 
point in 3D has no influence on a viscous fluid. 

By our physical common sense, all these results are easy to interpret, and we could guess to prove 
that a surface in 3D has a non- negligible effect on a viscous fluid as the curve in 2D, and that a curve in 
3D is an erasable singularity as the point in 2D for the Navier-Stokes equations. The goal of this article 
is exactly to show these two points. We will note that the first point can be obtained by a similar 
proof than [7J, but the case of the curve requires a different way to cut a divergence free function. 

More precisely, the motion of an incompressible viscous fiow in a domain 0, is governed by the 
Navier-Stokes equations: 

dtu- uAu + u-\7u = -Vp, \f{t, x) e {0,oo) X n, (1.1) 

where u = {ui,U2,U3) denotes the velocity, p the pressure and u the viscosity. The incompressibility 
condition reads 

div n = 0, V(t, x) G (0, oo) x Q. (1.2) 
The natural boundary condition is the no-slip boundary condition: 

u = 0, V(t,x) G (0,oo) X (90. (1.3) 

As we will study the asymptotic behavior when an obstacle C„ shrinks to a point, a curve or a surface, 
we need to give an initial data which allows us to compare solutions defined on different domains. For 
this reason, we use an important quantity for the 2D ideal flows (ideal means that = 0), which is the 
vorticity: 

uj := curln = {d2U3 - d3U2,d3Ui - diU3,diU2 - d2Ui). 

As in all the previous works |6l [71 [5] , we fix a smooth initial vorticity wq, which is divergence free 
and compactly supported in M^. Let {C^} be a family of smooth simply connected compact set of R^, 
then we claim that there exists a unique vector field Uq on r2„ := \ C„ such that: 

divuS = 0, curing} = wo I n„, t^o • ^lac„ = 0, u"^ £ L^{VLn). (1.4) 
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For such an initial velocity, it is well known that there exists a global weak solution n" of the Navier- 
Stokes equations p.ip - p.3|) on r2„ in the Leray sense (see Definition II .TP . and the main goal is to study 
the asymptotic behavior of u" when Cn shrinks to a point, a curve or a surface. 

In dimension two, the exterior of one obstacle is not simply connected, and the problem p.4p have 
an infinite number of solution in L^'°°. Indeed, we have to consider the harmonic vector field which 
behaves like l/\x\ at infinity. Then, to uniquely determine Uq, the authors in [SJ [7j have to fix the 
initial circulation: 

7 := ® Uq ■ T ds. 

JdCn 

In the case of one obstacle C„ (containing the origin) which shrinks homothetically to 0, the main 
result of [B] states that u"" converges to u in L2^^(M+ xM^), where u is the solution of the Navier-Stokes 
equations in with initial velocity uq verifying: 

divMo = 0, curl no = + 7^0) 

where 6q is the Dirac at the origin. Therefore, we do not feel the presence of the point, except in the 
initial data. In dimension three, the exterior of one obstacle is simply connected, and the circulation 
has no sense, so the reader should keep in mind the previous result in the case 7 = 0, i.e. a point in 
2D is erasable for the Navier-Stokes equations. 

When the obstacle shrinks to a Jordan curve T, [7j shows that the solutions converges to the 
unique solution in the exterior of the curve, where the initial velocity uq verifies 

divMo = 0, c\iT\uQ = ujQ + g^^^{s)6T, 

with (5r the Dirac on the curve T and Qi^^ a density which can be expressed in terms of wq, 7 and 
the shape of the curve. This function g^^^^ is continuous on F and blows up near the end-points of the 
curve. Actually, we can show that uq is a vector field continuous up to the curve, with different values 
on each sides, except near the end-points where it blows up like the inverse of the square root of the 
distance. Moreover, uq is tangent to the curve, and ^^jq is equal to the jump of the velocity through 
the curve. Then we see an important difference with the case of the point: even if initially 7 = 0, the 
function / which implies a non-negligible effect of the curve. 

A key point in dimension two is the explicit formula for the Biot-Savart law, law giving Uq in terms of 
Wo, 7, and the Riemann mapping T„, which sends to the exterior of the unit disk. In [6], the authors 
chose an homothetic convergence in order to write: r„(x) = ro(nx), and putting it in the Biot-Savart 
formula they obtained some good estimates. In [3, we assumed the convergence of r„ and T^, and we 
exhibited an example of where such a convergence holds. In that article, we were not be able to 
show that the convergence of the Riemann mappings holds for some "good geometrical convergence". 
The convergence of the biholomorphisms when the domains convergence is now established in [8j . 

To give the result in dimension three, we introduce the Biot-Savart law in M^: 

f ^ — y 

vq{x) = -1 — j ^xLJo(y)dy, (1-5) 

with X denotes the standart cross product of vectors in M^, then vq is the unique vector field in 
verifying: 

divtiQ = 0, curluo = cjq, vq G L^(M^). 
With this definition, we will prove that the solution Ug of ()1.4p appears to be the Leray projection (see 
later for the precise definition) of vq: 

K = IPn^^^o- (1.6) 
When one obstacle C„ shrinks to a point in dimension three, the authors in [5] shows that 

converges to a solution of the Navier-Stokes equations in with initial velocity vq. Therefore, a 3D 

viscous fluid does not feel a material point. 

The goal of this article is to treat the case of a curve and a surface in dimension three. 

To give a sense to a convergence of vector fields defined on different domains, we set the notation of 
extensions. For any function / defined on we denote by Ef the extension of / on M^, vanishing 
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on Cn- If / is regular enough and vanishes on dCn, then 'SJ[Ef) = E(Vf). Similarly, if is a vector 
field regular enough and tangent to the boundary, then div Ev = E(divv). 

The goal of Section [2] is to show a similar result to the case of the curve in dimension two. 

Theorem 1.1. Let C be a simply connected compact set o/M^, and ujq be a divergence free function in 
Cq°{M.^). We consider {Cn} a family of smooth simply connected compact sets converging to C (in the 
Hausdorff sense), such that C G Cn- Let u"" be a global weak solution to the Navier-Stokes equations 
on Qn '■= \Cn (in the sense of Definition \l. Tj j with initial velocity Uq ( (jl.6p the solution of (jl.4p ). 
then we can extract a subsequence such that Eu" converges to u in Lf^^(R+ X n), where u is a global 
weak solution of the Navier-Stokes equations on := \ C with initial velocity uq = fnvo- 

In this theorem, we do not need any assumption on the regularity of 17. We recall that the existence 
of a weak solution in non smooth domain is already known in literature (see e.g. \Tl). 

In the case where C is a smooth surface, we discuss in Subsection 12.31 about some properties of the 
initial datum, but it is clear from the no-slip condition that the surface is a non-negligible singularity. 

Section [3] is devoted to prove that the curve is an erasable singularity for the 3D Navier-Stokes 
equations. 

Theorem 1.2. Let T be a C'^ infective compact curve ofM?, and be a divergence free function in 
C^(^). We consider {Cn} a family of smooth simply connected compact sets converging to T (in the 
Hausdorff sense), such that F C C„. Let be a global weak solution to the Navier-Stokes equations on 
r2„ with initial velocity Uq (see (|1.6p ). then we can extract a subsequence such that Eu" converges to 
u m Lf^^(R+ xM^), where u is a global weak solution of the Navier-Stokes equations on with initial 
velocity vq. 

Therefore, we obtain exactly the same theorem than [S], but in the case of a curve, which is the 
similar result to the point in 2D (see [6]). If such a theorem seems natural, the proof requires a new 
way to cut divergence free test function. The method is based on the Bogovskii results. 

Notation and recall. 

We start by introducing the classical vector fields spaces for the study of incompressible flows. 
Definition 1.3. Let Q an open subset o/M^. We denote by 

• V{n) := |n G C^in) I divn = in l^}; 

• 'H{^) the closure ofV{Q) in the norm L^; 

• V{0,) the closure o/y(0) in the norm , and its dual space by V'{0,); 

• H2{n) := |u G L'^{Q) \ div u = in n, u ■ n = at 50 jQ ; 

• := jii G H^{n) I divu = in q|; 

• G{Q) := ju; G L2(0) I ti; = Vp, for some p G ///^^(O)}; 

Now, we give some results on these solenoidal spaces. The following theorem concerns the Leray 
projection [31 Theo III. 1.1]. 

Theorem 1.4. Let 17 be an arbitrary domain in M^. Then G{Q) and T-L{0,) are orthogonal subspaces 
in L^(ri). Moreover 

L'^iQ) = G{n)(B'H{n). 

Let us note that the previous theorem does not require any assumption of 0. The validity of such a 
decomposition in L*? requires some restrictions, whereas for L^, the Hilbert structure allows us to stay 
general. Moreover, we deduce from this theorem the existence of a unique projection operator (called 
Leray projection): 

¥n : L\n)^n{n), 



■ n should be understood in (see e.g. [3| Theo III. 2. 2]). 
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and by orthogonality we have: 

lliPf^-uIlL^ < \\u\\l^, yu G L^{n). 

Let us also mention that Theorem 11.41 implies the existence- uniqueness of a solution of (|1.4p , and 
that this solution is (11.611: 



Corollary 1.5. Let C be a simply connected compact set, and 0, := R \C. If ujq G ), then there 
exists a unique solution of 

curln = uJo\n in V'(0), u £ 
Moreover, u = P^^O) where vq is given in terms of ujq by (|1.5p . 



Proof. By Theorem II. 4[ we know that fo|n — Fn^o is a gradient and then curlP^fo = curluob = ^o\n- 
Then PqUq is a solution, which prove the existence part. 

Concerning the uniqueness, let us assume that u and w are two solutions. As = \ C is simply 
connected in R^, curl(n — w) = on Q implies that there exists p such that u — w = \7p. Then Vp 
belongs to Hq which is possible only for Vp = (see Theorem 1 1.4p . Its ends the proof. □ 

Actually, to infer that the previous corollary implies the existence- uniqueness of (|1.4p . we have to 
note that some spaces are similar. If T-l{Q) and H2{yi) seems to be the same (as V(0) compared to 
H^iyi)), we aware that there exist some counter examples, even in the case of smooth domains (see 
e.g. [3]). However, it will not be the case for the domains considered here. Indeed, in this article we 
always consider exterior domains 17 := \ C where C is a simply connected compact set. Then, [3l 
Rem III. 3. 4] implies that 

17 is a locally Lipschitz domain <^=^ Q. verify the cone condition. 

The reader can look for the precise definitions in [3J , but we state that the exterior of a smooth compact 
or of a thin surface (as in Subsection 12. 3p verifies the cone condition. 

Then, Theorems III. 2. 3 and III. 4. 2 from [3] can be applied to get the following. 

Theorem 1.6. Let Vt be the complementary of one simply connected compact set. If verifies the 
cone condition, then 

ni^) = H2{n) and v{n) = H\n). 

The exterior of the curve verifies also the cone condition, but u ■ n has no sense in this case. 
Putting together Corollarv 11.51 and Theorem we finally get, in the case of smooth simply con- 
nected compact obstacle Cn, the existence- uniqueness of a solution of ()1.4p . and that this solution is 

(USD. 



Now we give precise definition of a global weak solution of the Navier-Stokes equations in the Leray 
sense. 

Definition 1.7. Let uq G 7^(17). We say that u is a global weak solution of the Navier-Stokes equations 
on 17 with initial velocity uq iff 

• u belongs to 

c{[o, oo); v'm) n L^,{[o, ooy,nm) n lI,{[o, oo); Vim 

• u verifies the momentum equation in the sense o/V'(17), i.e. \/^p G C"'^([0, oo); V(17)), we have 
for all t: 

[ {u-ij){t,x)dx + [ [ {-u-il^t + vVu -.Vi) - [u^u) ■.Vi)){t' ,x)dxdt' = [ uo-i^iO,-); (1.7) 
Jn Jo Jn Jn 

• u verifies the energy inequality: 

i2(Q) < ||'Uo|li2(!^) Vt > 0. (1.8) 
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Without any assumption about the regularity of 17, the Leray theorem states that there exists a 
global weak solution of the Navier-Stokes equations in the sense of the previous definition (see e.g. [U 
Theo 2.3] and Theo III.3.1]). So the goal here is to choose one solution it" on r2„ and to prove that 
the limit of is a weak solution of the Navier-Stokes equations on 17 or on M^. 

2. Convergence of viscous flows when the obstacles converge to a singular compact 

2.1. Convergence of the initial velocity. 

Let C be a simply connected compact set. We consider {C„,} a family of compact, simply connected, 
smooth set of such that C C Cn and Cn converges to C in the Hausdorff topology when n — )■ oo. 
The only property needed here concerning Hausdorff topology is the following. 

Proposition 2.1. // {Cn)ne'N converges in the Hausdorjf sense to C and K is a compact subset of 
O := M'^ \ C, then there exists uq such that Qn := M'^ \ C„ D -fC for all n > uq. 

As mentioned in the introduction, we fix LOO G Vi^^) and we define Uq as in (jl.6p . which is the 
unique vector field verifying (jl.4p . Let us show that Uq converges strongly to uq := WqVq. We recall 
that E is the extension by zero in C„. 

Proposition 2.2. With the above notations, we have that 

Euq — )• no strongly in L'^{Q). 

Proof. As the Leray projection is orthogonal in L^, we get from (jl.6p that 

ll^'^o = ||'u-oIIl2(q^) < ||wo||l2(!^,^) < ll-yollLacn)- 

By the Banach-Alaoglu's theorem, we infer that there exists wq G and a subsequence n — )• oo, 

such that 

Euq Wq weak in L'^{Q). 
This weak convergence implies in particular that 

II'^^'oIIl^ ^ liiiiinf llE'uo 11^2, 

and 

cmlvjo = ujqIq in V'(il). 

In the proof of Corollarv 11.51 we prove that curlPQ^t^o = curl7Do|r2„ = ^olfin- Then, Corollarv 11.51 
implies that Uq = Fq^wq and then 

||£^?iollL2 < II-woIIlz. 

Putting together the two last inequalities, we get ||SnQ||j;^2 — )• ||ii;o||L2. Using the weak convergence in 
L^(r2) of Euq to f^o we obtain the strong convergence in L^(0). 

As Eul^ belongs to n{n n), we deduce directly from the strong convergence that wq belongs to 
T-L^Q). Then we have two functions in T-L{yi) having the same vorticity, which implies that uq = wq (see 
Corollary II. 5p . 

□ 

Remark 2.3. We have obtained in the previous proof that 

curl no = Wo = curlwo in V'{VL) 

but we do not have 

curl no = Wo = curl no in H^'^i^"). 

Even in the case where S" is a thin surface, we can just pretend that curl (no — no) belongs to H~^{F?) 
and supported on the surface. Actually, we will prove in Subsection 12.31 that 

curl no = iOQ + gs5s 

where gs is the jump of the tangential component of no. 

In the case of the curve, we will show in Subsection 13.11 that there is no function belonging in 
if~^(M'^) compactly supported on a curve, and we will obtain that no = no- 
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2.2. Time evolution. 

For all n, we denote by u" a global weak solution of the Navier-Stokes equations (in the sense of 
Definition II. 7p on fi^ with initial data Mq. 

By Proposition [221 we already know that Euq — )• uq in Moreover, thanks to the energy 

inequality (jl.Sp . we state that 

Eu" is uniformly bounded in L°°((0, oo); "Hq) and VEvT" is uniformly bounded in L^((0,cxd) x ^^). 

(2.1) 

Now, we need to establish a temporal estimates. First, we use that verifies the Dirichlet boundary 
condition for t > in order to write the following Sobolev inequalities: 

< C\\Eu-\\%,^\\V{Eun\\%\^s) = C\^^^^^^^ (2.2) 

where C is independent of n. We consider in a first step O an open bounded set relatively compact 
in Q. By Proposition 12.11 there exists no such that O f^Cn = ^ for all n > no, then we have for any 
$ G V{0) and n > no'- 

|(u"(t),$) - (^"(s),^)! = [ [ -i^Vn": V$ + («"«) n") :V$ 

S) \\U |lioo(K+,L2(Q^))|| VU. lli2(K+xn„)ll^^ll-f^2 

< C((t-.)i/2 + (t-.)V4)||$||^,, 

where we have used (j2.2p . This inequality implies that is equicontinous as a family of functions 

from M+ to V'(0). Using that this family is bounded in L°°(M+; L^(0)), Ascoli theorem gives that 
{Eu''} is precompact in L'^{R+; H-'^{0)). Moreover, {^u"} is also bounded in Lf^^{R+;H^{0)), 
then we get by interpolation that this family is precompact in L^q^(M+ x O). Finally, by a diagonal 
extraction on the compact sets of we conclude that 

{£'u"} is precompact in Lf^^{M.'^ x Jl). 

Therefore, there exists a subsequence n — )• cxd such that Eu^ — )• u strongly in Lf^^(M.'^ x Q). Moreover, 
extracting a subsequence if necessary, we know from ()2.ip that 

u G L°°((0, oo);n{n)) and Vn G L2((0, oo) x Q). (2.3) 

For any test function ip G Ci(M+; V(l^)), there exist V''' G C^{{0,oo) x J]), div^^ = such thalS 

V'*^ ^ V strongly in L?„,(M+; V(J^)), Vg G (1, oo). 

For k fixed, there exists such that supp ^/^^ n = 0, Vn > rijt, then ()1.7p reads 

(-u" • V'f + i^Vu" : VV''' - (n" u") : VV''') dx dt' = 0. 

Thanks to the strong convergence of in L^^^(R^" x VL) and the weak of Vn" in L^(M^" x J7) then we 
can pass to the limit n — )• oo to get: 

{-u ■ V'f + vVu : VV'^ - (u u) : VV''') da; dt' = 0. 

i 

In particular, this equality putting together with ()2.3p and ()2.2p gives that 

'4/3 
^loc 



belongs to Lf/^^(]R+, V'(Sl)). 



Then we can pass to the limit in 

(dtu ■ + uVu : V^'' -{u0u): V^)^) dx dt' = (2.4) 



see e.g. [3 Prop 3.6] 
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as — )• OO to get 

Mif G C^(0, oo), / / {dtu ■ ipip + uipVu : Vip - ip{u u) : Vip) dx dt' 
Jo Jn 



0. 



This equality implies that 



^ / u- ip = I {u-ipt- vSIu -.Vip + {u(E>u) : Vip) dx (2.5) 



n Jn 



in the sense of distribution in M+. Since the right hand side term belongs to Lj^^(M+), the equality 
holds in Ll^(R+). 

Now, we check that u belongs in the good functional space. Thanks to the previous equality, and 
, then we can easily prove that 



u E c{[o,oo);V\n))nc^{[o,oo);n{n)). 

This argument can be found in [9, Subsection III. 3.1]: as dtu belongs to -L^(V'), then its implies that 
u is almost everywhere equal to a function continuous from into V'. Moreover, using the fact that 
u G L°°('H), then [9l Lem 1.4] states that the continuity in V' implies the weak continuity in time with 
values in H. 

Moreover, thanks to the continuity in V', we infer that the equality ()2.5p in the sense of Lj'Q^(M^) 
implies that the integral equality ()1.7p holds for all t > 0. Indeed, for the initial data we know from 
the uniform convergence in iifj~^(il) that Euq — )• u\t=o in H^'^{Q). However, we proved in Proposition 
12.21 that Euq —7- uq in L^(r2), which allows us to state by the uniqueness of the limit in that the 
initial velocity is uq = Pf^^^o- 

To finish the proof of Theorem II. 1|, we have to prove the energy inequality. We take the liminf of 

limjnf||i?n"(t)||2,(^)+2i.liminf^ \\VEu^{t)\\l,^^^ < hoWl^^^^ 

and we note that the weak limii0 in of Eu^{t) to u{t) and the weak hmit in -L^((0, t) x of VEu"' 
to Vn imply that 

lhWlli2(n) <liminf ||£;u"(t)||2 2(^^ and || VuII^^^^q^^^^^^ < liminf HV-En^H^^^^Q^^^^^^. 
It gives the last point required in Definition 11.71 which ends the proof of Theorem 11.11 

Therefore, we have shown that the Navier-Stokes solutions converge when domains convergence to 
the exterior of a simply connected compact set C. We note here that we do not assume any assumption 
on the regularity of C. In particular, this result holds if C is a surface S. As u{t, •) G V(ri) for a.e. time, 
it is clear that the surface has a non-negligible effect on the motion of 3D viscous flow: u verifies for 
almost every time the no slip boundary condition. In the following subsection, we discuss about the 
initial velocity properties in the particular case of the surface, and the goal is to get some similarities 
to the curve in 2D. 

2.3. Remark on the behavior of the initial velocity near a smooth surface. 

In the two dimensional case, we obtain in [7] an explicit formula of the initial velocity in terms of the 
biholomorphisms. This formula allows us to state that no is continuous up to the curve with different 
values on each side, except near the end-points where it behaves as the inverse of the square root of 
the distance. In our case, we do not have this formula, and we use in this subsection classical elliptic 
theory in order to get similar results in the case where 5 is a bounded orientable surface of codimension 
1 in R^. 



"^Indeed, by the uniform estimates and diagonal extraction, we can find a common subsequence such that Eu"{t) 
weakly converge for all t G Q+. Then, we conclude by the continuity that this sequence holds for all t G R . 
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As mentioned in Remark 12. 3|, uq — vq is curl free in il. and as is simply connected we infer that 
there exists p such that uq — vq = Vp. We know that vq is continuous up to the boundary, then the 
goal is to determine the behavior of Vp near S where p verifies the following elliptic problem 

Ap = inn 

dp 

— = -vq -n on S, 
on 

where vq is regular on S. 

This subsection is independent of the convergence theory, and the goal here is just to give an 
example of behavior of uq. Therefore, we add here some conditions on S in order to apply classical 
elliptic results. We assume that S* is a C°° manifold and its boundary ;B is a closed curve. 

The study of elliptic equations in the exterior of a surface with the Neumann condition is standard 
for the crack problem in 3D linear elasticity. Actually, to get exactly the Neumann boundary condition, 
we add a regular function h such that p : = p + h verifies 

dp 

Ap = f in n, — = on 5. 

on 

For ojQ regular enough, we rich the necessary regularity for / in order to state that p have an expansion 
near B on the form 



k>0 



in local polar coordinates {r,9). Such a result is proved in [2j (see also the references therein). In 
particular, it implies that uq is continuous up to S, with possibly different values on each side, except 
near the boundary B where uq behaves like the inverse of the square root of the distance. Therefore, 
we obtain exactly the same behavior in 3D in the exterior of a surface than in 2D in the exterior of a 
curve. 

Thanks to the continuity up to the surface, it is easy to see that the tangent condition implies that 

div no = and curl uq = loq + gds 
in V'(M^), where g is the jump of the tangential component of uq through the surface. 



3. Viscous flow around a curve 

As in the previous section and as in [5J, we consider C.„ a family of compact, simply connected, 
smooth set of such that Cn converges to T in the Hausdorff topology when n — )• oo. Here, T is 
assumed to be a compact injective curve, i.e. of dimension 1 in M^. The goal of this section is to 
prove that the curve is a negligible singularity for the Naviers-Stokes solutions in M^. 

3.1. Convergence of the initial velocity. 

In the following subsections, we will need to use a suitable cutoff function of a small neighborhood 
of r. Let X be a function verifying 

XGC°°(M), x(s) = Oon(-oo,l), = 1 on (2, +oo), 

then we define 

r]n{x) := x{nd{x,T)y (3.1) 

Although it is obvious that r/„ vanishes in a small neighborhood of F, we need the following estimates 
of ||Vr/„||iP. 

Proposition 3.1. Let T be a C'^ compact injective curve. There exist no and C > such that for all 
n > no, 

l|V^n||Loo(R3) < Cn and mes^supp (1 — ??n)^< C/n^, 
where mes is the Lehesgue measure. 
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Even if such a proposition seems standard, we write the details in Appendix [A] for a sake of com- 
pleteness. Thanks to this proposition, we can prove that there is no function supported on T which 
belongs to H~^{ 



Lemma 3.2. Let f a function belonging in H ^(M^). If 

{f,^)H-i,H^=o, v^Gcr(M3\r) 

then f = in V'{R^). 

Proof. We fix (y9 G C^(M^), and we introduce 

where r/„ is the cutoff function defined in ()3.ip . As we proved in Appendix El that rjn is supported in 
the exterior of a small neighborhood of T, then by assumption (and density of C^(R'^ \r) in C^(M^ \r) 
for H^{W^) norm), we have that 

{f\^n)H-\m = 0, Vn. 
Thanks to Proposition I3.H we get easily that 

\\ipn - ^p\\l2 < C{ip)/n and \\V(pn - Vv?||2,2 < C{(p). 

By Banach-Alaoglu theorem, we can extract a subsequence such that (/?„, — ^ weakly in H^(U.^). 
In particular, it implies that 

(/, V?) = if, Vn) + Vn) = {f,V- fn) ^ 

as n — )■ DO. Then, we have proved that for all 99 S C^(M^), we have (/, ip) = 0. □ 

Now, we can come back to our main problem. Let T be an injective compact smooth curve in M^. 
We consider C„ a family of compact, simply connected, smooth set of M.^ such that Cn converges to T 
in the Hausdorff topology when n — )• 00. Let loq an initial vorticity which is bounded and compactly 
supported in W^. As in the previous section we set 

:= i^ob,, with Qn ■= \ Cn 

then we define by Uq := Pn„vo the unique vector field in r2„, solving ()1.4p . We prove in the following 
proposition that we do not feel the presence of the curve for the limit initial velocity. 

Proposition 3.3. With the above notation, we have 

Euq vo strongly in L^(M^), 

where vq is the velocity field without obstacle ()1.5p . 

Proof. As vq is continuous and behaves like 0(1/ at infinity, we obtain directly that vq belongs to 
LP(M3) for any p E (3/2, 00]. 

Firstly, we introduce the stream function corresponding to vq: 

f X — y 

^(x) = - / — X voiy) dy, 

Jr3 47r|x - y]-^ 

which verifies curl?/; = vq, divtp = and HV'IIl"" < C'||?;o|| 2,2^2,4 (see e.g. fS]), so ^p is bounded. 
Secondly, we define 

w"" := curl (ijnijj) = rinVQ + Vr?„ x ip, 
where r/„ is the cutoff function (|3.ip . By construction, we deduce that w'^ belongs to H2{^n) = ^(^n)- 
Finally, we use that Uq is the projection of vob„ on T~{-{0,n), to compute 

ll^'"0 - ^o||l2(ib;3) < \\uq -Vo\\L2(n„) + \\vo\\L2(c„) 

< Ww"" -vo\\L2(^nn) + \\vo\\LHCn) 

< \\Vr]n X i'h^^n^) + 11(1 - f?n)^'o||L2(f^^) + ||?;o||l2(c„) 

< l|Vr/„||i2||?/;||Loo + ||t'o||L°°(l/n + mes(Cn)) 

< c, 



10 



C. LAC AVE 



where we have used Proposition 13.11 This inequahty ahows us to pretend that there exist v £ L^(M^) 
and a subsequence such that 

Euq — vo ^ V weakly in L^. 

By passing to the hmit, we obtain that divv = curlu = for all test function in C^(M^\r). Moreover, 
divf and curlt; belongs to H~^{^'^) and we apply the previous lemma to state that divi; = curlu = 
on M^, and then v = 0. 

Moreover, reasoning as in the proof of Proposition \2.2\ we pass from the weak convergence to the 
strong convergence of Uq to vq, which ends this proof. □ 

3.2. Proof of Theorem 

The begin of the proof follows the same idea of the proof of Theorem ll.lt we prove that we have a 
strong limit of £'n" to u a solution of the Navier-Stokes equations in "the exterior of the curve" . The 
second step is to show that u is actually a solution in M^. 

For all n, we consider u" a global weak solution, in the sense of Definition 11.71 of ()l.ip - ()1.3p in 
with initial datum Mq, which is defined in the previous subsection. 

By Proposition 13.31 we already know that Euq — )• vq in Next, we do exactly the same thing 

than in Subsection 12.21 

• uniform estimates of ^n" in L°°{{0, oo); L'^{R^)) and of VSu" in ^^((O, oo) x R^); 

• equicontinuity in V^^(O), with O relatively compact in \ T; 

• precompactness of {Eu^} in Lf^^(R'^ x (M^ \ P)^; 

• by a priori estimates we know that 

u G L°°((0,oo);L2(M3)) and Vn G L2((0,oo) x R^); 

• passing at the limit n — )• oo in the weak formulation of the momentum equation: V^/^*^ G 
C^((0,oo) X (M3 \ r)) I div V'' = 0, we have 

Now, we need to prove that the momentum equation is verified for all ^jJ G C^(M'''; V(M'^)). Then, 
to finish the proof of Theorem 11.21 we have to establish ()1.7p for test functions which meet the curve. 
The following lemma will be the key of this extension. 

Lemma 3.4. For all if G C^((0,oo) x M'^) such that div if = 0, then there exists a sequence tpn G 
C^°°((0,oo) X (M^ \r)) such that 

div^/;„ = andipn^V weak-* in L°°(M+; 

Proof. All the difficulty comes from the condition divipn = 0. Indeed, without this condition, it is 
sufficient to multiply by the cutoff function. In [6l[7l[5], the standard way to construct divergence free 
functions compactly supported outside the obstacle is to multiply the stream function by the cutoff 
function. However, we see in the proof of Proposition 13.31 that the computation of the norm makes 
appear ||V^?7n||i,2 which blows up strongly in our case. 

So, we present here a new way to approximate divergence free function, which is not explicit as in 
the standard way. The following method comes from [3l, Chap III. 4] and is based on the Bogovskii 
operator. We set 

fn{t,x) := r]n{x)ip{t,x) 
which belongs to C^((0, oo); C^(M^ \ P)) but which is not divergence free. Its divergence 

gn{t,x) := d\\ {rin{x)Lp{t,x)) = ip(t, x) ■ Vr]n{x) 
verifies the following estimates for all t: 

/ 9nit, •) = and\\gn{t, ■)\\lp{r3) < Cn^p, 

where C depends only on 93. To correct this divergence, we use the result of Bogovskii. Let B he a. 
ball big enough containing the support of As the set ^ := B\T verifies the cone condition, and as 
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gn G -L^(^i) with /^^n = 0, then Theorem 3.1 in [3l Chap III. 3] states that there exists at least one 
solution hn of the following problem: 

divhn = -gn, /l„ G Wo'^(r2), WhnWw-^-P < Cp\\gn\\LP- 

Moreover, the constant Cp depends only on il. and p, and as gn as a compact support in il., so is 
Then, we can define for all t, g-nit, •), and we have that hn ^ strongly in L°°(M+; W^'P{M.^)) for all 
p < 2 and uniformly bounded in ; H (Mr)) . S o we can extract a subsequence such that hn 

converges weak-* in L°°(R^; //^(M'^)) and by uniqueness of the limit in L°°(M+; W^'^{B)) it converges 
to 0. 

From the definition we get that fn — f converges to strongly in L°°(M+; ^^^'^(M^)) for any p < 2 
and weak-* in L~(M+; ^^(M^)). 

In conclusion, ipn ■= fn + hn is divergence free, compactly supported in i? \ F and converges to if 
weak-* in //^(M^)). Smoothing ipn by some mollifiers, we obtain the result. □ 

In the previous proof, we see that the approximation constructed verifies also 

dtiJn dtipweak-* in L°°(M+; ^^(M^)). 

Therefore, we apply the previous lemma in order to pass to the limit in (j2.4p . which implies that this 
equality is verified for all test function in C^((0, oo) x M^). Now, we finish as we did in Subsection 12.21 

• (EaD gives that dtu G Lfl^(R+;V'{R^)); 

• we pass at the limit to say that the momentum equation is verified for all ip G C^(M"^; V(M^)) 
in LUR^); 

• we get the continuity in time with values in V'(M'^); 

• we conclude to the validity of ()1.7p for all t; 

• identification of the initial velocity: n(0, •) = vq by Proposition 13. 3| 

• thanks to weak convergence, we prove the energy inequality (jl.Sp on M^. 

Its ends to prove that n is a global weak solution of the Navier-Stokes solutions in M^, in the sense 
of Definition 11.71 

We finish the proof of Theorem 11.21 by noting that the uniform estimates ()2.ip and the convergence 
in L'^ocO^'^ ^ (^^ \ r)) implies that 

Sii" u strongly in Lfo^(M+ x M^). 

Indeed, for any e > and K a compact of M"*" x M^, we decompose the norm thanks to ()2.2p : 

WEu'' - u\\l^k) < WVkiEu'' -u)\\L2^K) + \\i'i^-m)iEu'' -u)\\l2^k) 

< \\Eu" - u||i2(^nsupp + 11(1 - m)\\L4Eu'^ - u\\]^^^^2)\\Eu'^ - u\\%\o,T;m)'^^^^ ^ 

where T is a time bigger than all t €z K. Using the uniform estimates, we deduce from Proposition 13.11 
that the second right hand side term tends to zero as A; — )• oo, uniformly in n. So there exists ke such 
that this term is less than e/2. For this kg fixed, the first right hand side term tends to zero as n — )• oo 
because the compact does not meet F. Then, there exists such that this term is less than e/2 for 
any n> Nf,, which ends the proof. 

Appendix A. Cutoff functions 

In order to differentiate the function x i— )• d(x,F), we have to check that the minimal distance 
d(x, y) is reached for a unique y G F, at least for x closed enough to F. We need to assume that the 
curve is at least C^, because such a property is false near a corner. For example, if F is the curve 
{(x,y,0)|x G [—1, l],y = |x|} then for all M on the half line {(0, y, 0)|y > 0}, we have two minimums 
on F: Ml := {-y / V2,y / ^/2,{)) and M2 := (y/\/2, y/^/2, 0) and d(M,F) = d(M,Mi) = d(M,M2) = 

y/V2. 

This example shows that we have to assume some regularity for the curve. 
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Lemma A.l. Let T be a injective curve, then there exists no such that 

Vx I d(x,r) < 1/no, 3!y e r I d{x,y) = d{x,T). 

Proof. Let s i— )■ r(s) = M he the arclength parametrization. Then, we recaU the standard definition: 

• the tangent vector and the curvature 

^(s) = ^ = p{s) = = l|r"(s)||; 

• the main normal vector and the binomial vector (if p[s) 7^ 0) 

As r is assumed to be C^, there exists pM > such that \p{s)\ < pM on T. 

The existence of y G F such that d(x, y) = d(x, F) is obvious because the map z 1— )• d{x, z) is 
continuous on the compact F. Let us assume that the conclusion of the lemma is false, i.e. that the 
infimum is reached twice: 

Mn £ N, 3xn I d(x„,F) < l/n, 3y„ = F(si,„),z„ = F(s2,„) 

such that 

si,n < S2,n and d(x„,y„) = d(x„,F) = d(x„,2;„). 

Extracting a subsequence, we have by compactness that Xn,yn,Zn — )■ x G F. As the curve is assumed 
to be injective (so without cross-point), we infer that — S2,n — ^ 0. So, by continuity of F in x, we 
note easily that there exists N such that 

d(xiv,F(s)) < l/(2pA/), Vs e [si,jv,S2,iv]. (A.l) 
Next, we introduce the following function 

f:s^ d{xN,r{s)f = \\xNM{s)f, 

and we differentiate twice to get 



f'{s) = 2xnM{s) ■ ^{s), f "{s) = 2(1 + p{s)xnM{s)-1){s 



The previous computation holds even if p{s) = 0. 
By assumption, / is minimal for si^n and S2,Af) so 

/'(S1,JV) = = /'(S2,JV), 

but (jA.ip implies that 

f"{s) > 1, VS e [si,jv,S2,iv] 

which is impossible. This contradiction allows us to end this proof. □ 



Remark A. 2. In the previous proof, we infer from /'(si^at) = that xnM{si^m) ■ '^{si,n) = when 
M{si^n) is the closest point of F to xn- 

Now, we link the differentiability of g : x 1— )■ d(x, F) and the previous lemma. 

Lemma A. 3. Let T be a C"^ injective curve and fiQ from Lemma \A.1[ Then for all x £ M.'^\T such 
that d(x,F) < 1/no; 9 is differentiable in x and 

Vg{x) 



If - y\\ 

where y is the unique point ofV such that \\x — y\\ = d(x,F). 



3D VISCOUS INCOMPRESSIBLE FLUID AROUND ONE THIN OBSTACLE 13 

Proof. This proof can be found in the exercise book [4J and we copy it for a sake of completeness. 
We fix X G \ r such that d(x,r) < 1/no, then Lemma [A.ll states that there exists a unique y £ T 
verifying ||x — y|| = d(x,r). 

The first step consists to show the fohowing property: for all /i G M'^ small enough, we set yh by the 
unique point of T such that Ijx + /i — = d(x + h,T), then 

- y|| as /i 0. (A. 2) 

Let us assume that (jA.2p is false, then there exist a sequence hn ^ and 6 > such that \\yhn —y\\ > S 
for all n. We introduce the following compact set 

G:={zer\ \\z-y\\ > 6} 

and the function 

h : z € G \\x — 

This function is continuous and reaches his minimum. As y ^ G, the uniqueness part of Lemma lA.ll 
allows us to check that this minimum is strictly greater than ||x — Hence, we have 

d{x,G) > d(x,r). 

However, we use the fact that y^^ belongs to G to write 

g{x + hn) = d(x + hn, T) = \\x + hn - yh-a II = d(x + G) 

implying that 

g{x + hn) d{x, G) > g{x) 
which is impossible by the continuity of g. The property ()A.2p follows from this contradiction. 

Let us prove in the second step that (|A.2|1 implies that g is differentiable. On the one hand, we have 
that 

g^{x + h)<\\x + h- yf = g'^{x) + 2{x-y)-h+ \\hf. (A.3) 
On the other hand, we compute 

g^{x + h) = \\x + h-yhf = \\x-yhf + 2{x-yh)-h+\\hf 

> \\x - yf + 2{x-y)-h + 2{y - yn) ■ h + \\hf. 

Hence, ()A.2p implies that 

g'^{x + h)> c/2(x) + 2{x-y)-h + o{h). (A.4) 
Putting together (|A.4p - (|A.3p . we conclude that g"^ is differentiable and that 

V(5'(x)) = 2g{x)Vg{x) = 2{x - y) 
which ends the proof. □ 

From the value of Vg, we see that the uniqueness part of Lemma lA.ll is mandatory to obtain the 
differentiability of g (or of g"^). 

Finally, we can establish the estimates of the cutoff function by proving Proposition 13.11 



Thanks to the definition and Lemma lA. 31 we compute for n > 2no 

V?7„(x) = nx'^nd(x, r)^ Vd(x, r) = nx'^nd(x,r)^ ^ 



\x - y\\ 

where y is the unique point of F such that ||x — y\\ = d(x,F). Indeed, x' ('^d(x, F)^ / iff 1/n < 

d(x,F) < 2/n, so X G \ F and d(x,F) < 1/no which allows us to apply the previous lemmas. 
Hence, we obtain directly the first point: 

l|V^ri||L°=(M3) < "-||x'IIl°°{R)- 

Concerning the support, we use Remark lA. 21 to note that for all x G supp (1 — r/„), then the unique 
y G F such that d(x, y) = d(x, F) verifies • = 0. If there is an interval where p{s) = 0, then it 
implies that F on this part is a segment. Then, it is obvious that the volume such that d(x, F) < 2/n is 
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©(l/n^) on this section. So, without loss of generahty, we assume that p{s) / on F, and the previous 
remark imphes that 

supp {l-r]n) = {T{s) + r cos (flJ'is) + r sin 9?^(s) \{s,r,ip) e [0, L] x [0, 2/n] x [0, 27r)} |J O, 

1=1,2 

and the appUcation ^: {s,r,ip) i— t- r(s) + rcos(/?7^ + r sin 99"^ is a diffeomorphism invertible where 
^~^(M) gives the parameter of T{s) where d(M, F) is reached, the distance d(M, F), and the angle 

between T{s)M and 1^. In the expression of supp (1 — rjn), Oi denotes the half ball on each endpoints, 
which obviously verify mes Oi = 0{l/n'^). 
Using the general relation 

^ = e{s)i^{s), ^ = -p{s)i'{s) - e{s)^{s), 

as as 

where is a smooth function called the torsion, then we can compute in the orthonormal base 

(I — rp cos ip 
r9 sin if cos if —r sin ip 
—rO cos <p sin Lp r cos 99 

which implies that | det L'^'(r, s, /)| < r/2 if r is chosen small enough (i.e. r < l/(2pAf)). Finally, we 
can compute for n big enough: 

mes (supp (1 - r/„)) = / Idxdydz = f f f \det D'i?{r,s, f)\dipdrds + O {I /n^^ 

< / r/2dipdrds + 0{l/n^) = ^ + 0{l/n^) 

Jo Jo Jo " 

which ends the proof of Proposition 13.11 

Acknowledgements. I am partially supported by the Agence Nationale de la Recherche, Project 
MathOcean, grant ANR-08-BLAN-0301-01 and by the Project "Instabilities in Hydrodynamics" fi- 
nanced by Paris city hall (program "Emergences") and the Fondation Sciences Mathematiques de 
Paris. 



References 

[1] Chemin J-Y., Desjardins B., Gallagher I. and Grenier E., Mathematical Geophysics: An introduction to rotating 

fluids and to the Navier-Stokes equations, Oxford University Press, 2006. 
[2] Costabel M., Dauge M. and Duduchava R., Asymptotics without logarithmic terms for crack problems, Comm. 

Partial Differential Equations 28 (2003), no. 5-6, 869926. 
[3] Galdi G. P., An introduction to the mathematical theory of the Navier-Stokes equations. Steady-state problems. Second 

edition. Springer Monographs in Mathematics. Springer, New York, 2011. 
[4] Gonnord S. and Tosel N, Calcul differentiel. Ellipse, 1998. 

[5] Iftimie D. and Kelliher J., Remarks on the vanishing obstacle limit for a 3D viscous incompressible fluid, Proc. Amer. 
Math. Soc. 137 (2009), no. 2, 685-694. 

[6] Iftimie D., Lopes Filho M.C. and Nussenzveig Lopes H.J., Two Dimensional Incompressible Viscous Flow Around a 
Small Obstacle, Math. Annalen. 336 (2006), 449-489. 

[7] Lacave C., Two Dimensional Incompressible Viscous Flow Around a Thin Obstacle Tending to a Curve, Proc. Roy. 
Soc. Edinburgh Sect. A, Vol. 139 (2009), No. 6, pp. 1138-1163. 

[8] Lacave C., Two Dimensional Incompressible Ideal Flow Around a Small Curve, to appear in Comm. Partial Differ- 
ential Equations. 

[9] Temam R., Navier-Stokes Equations, Theory and Numerical Analysis, North-Holland, Amsterdam, 1979. 

(C. Lacave) Universite Paris-Diderot (Paris 7), Institut de Mathematiques de Jussieu, UMR 7586 - CNRS, 
175 RUE Du Chevaleret, 75013 Paris, France 
E-mail address: lacave@math.jussieu.fr 



